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We have used two-fluid dynamics to study the discrepancy 
between the work of Thouless, Ao and Niu (TAN) and that 
of lordanskii. In TAN no transverse force on a vortex due to 
normal fluid flow was found, whereas the earlier work found 
a transverse force proportional to normal fluid velocity u n 
and normal fluid density p n . We have linearized the time- 
independent two-fluid equations about the exact solution for 
a vortex, and find three solutions which are important in the 
region far from the vortex. Uniform superfluid flow gives rise 
to the usual superfluid Magnus force. Uniform normal fluid 
flow gives rise to no forces in the linear region, but does not 
satisfy reasonable boundary conditions at short distances. A 
logarithmically increasing normal fluid flow gives a viscous 
force. As in classical hydrodynamics, and as in the early 
work of Hall and Vinen, this logarithmic increase must be cut 
off by nonlinear effects at large distances; this gives a viscous 
force proportional to u n /lnu n , and a transverse contribution 
which goes like u n /(ln it n ) 2 , even in the absence of an ex- 
plicit lordanskii force. In the limit u n — * the TAN result 
is obtained, but at nonzero u n there are important correc- 
tions that were not found in TAN. We argue that the Magnus 
force in a superfluid at nonzero temperature is an example of 
a topological relation for which finite-size corrections may be 
large. 

PACS: 67.40.Vs, 47.37.+q 



I. INTRODUCTION 

In recent years papers have been published which ap- 
pear to give an exact expression for the transverse force, 
including the superfluid Magnus force, acting on a quan- 
tized vortex moving in a neutral superfluid These 
methods were extended to give the transverse force on a 
flux line in a superconductor in the ultraclean limit 
and to give a corresponding (but not easily calculable) ex- 
pression for the longitudinal force on a vortex in a neutral 
superfluid jl . These results have been criticized on 
the grounds that they do not agree with the widely, but 
by no means universally, accepted results in the litera- 
ture ||||. In particular, a transverse force proportional 
to the superfluid density, and independent of the nor- 
mal fluid density and velocity, was found, whereas it is 
generally supposed that the asymmetrical scattering of 
phonons or rotons by a vortex, demonstrated many years 



ago by Lifshitz and Pitaevskii jHJ] and by lordanskii jlTj , 
should lead to a contribution to the transverse force pro- 
portional to the normal fluid density and to the velocity 
of the vortex relative to the normal fluid component. 

In order to shed new light on this controversy we have 
decided to examine carefully the predictions of the two- 
fluid model for a pinned vortex in a moving fluid. The 
two-fluid equations should be valid under conditions in 
which the pressure, temperature, and normal and super- 
fluid velocities are varying slowly on the length scale set 
by the normal fluid mean free path. This cannot be true 
close to the core of a vortex, but it can be true every- 
where if the vortex is pinned by a sufficiently large cylin- 
der, on which either zero velocity or zero rate of shear 
(slip) boundary conditions are imposed. The equations 
are also valid at sufficient distances from the pinning cen- 
ter, even if the pinning center has a radius less than the 
mean free path. This line of thought is not entirely new, 
since similar arguments were used in the original work of 
Hall and Vinen fl^| on vortex arrays in helium, and by 
Vinen Jl3| on mutual friction. Their arguments, however, 
did not address the question of a transverse force. 

In sec. [n] we write down the equations of two-fluid the- 
ory for steady state flow, and examine the simple equi- 
librium solution associated with a vortex in a stationary 
fluid. 



In sec. Ill we set up linearized steady state equations 
and study their solutions. We want to find solutions lin- 
earized about the equilibrium vortex solution, and ex- 
amine the behavior of the steady state solutions at large 
distances from the vortex. Such an approach should be 
adequate to derive any transverse force proportional to 
normal fluid velocity and superfluid circulation, such as 
is quoted in refs. |[ |^ and |?], as well as to derive the force 
proportional to superfluid velocity and superfluid circu- 
lation. The partial differential equations are separable in 
cylindrical polar coordinates, so by a combination of an- 
alytical and numerical methods we can get a fairly com- 
plete description of those solutions that have the same 
symmetry as uniform flow past the vortex. 

In sec. IV we examine the influence of boundary con- 
ditions on the solutions of the equations of motion. In 
our initial discussion we look at the case of a cylindri- 
cal inner boundary on which the normal fluid is at rest 
and the temperature is uniform. There is a complica- 
tion, known long ago for the classical theory of flow of 
viscous fluid past a cylindrical obstacle, and discussed 
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for superfluids in the works of Hall and Vinen [12] and 
Vinen 13 1, which is that the linearized equations are not 
valid at very large distances, where the Reynolds number 
is of order unity even though the speed of flow is small 
Jl4| . This is true for any reasonable set of short-distance 
boundary conditions. This difficulty was circumvented 
by Oseen, and we can use the same technique. The drag 
force, in the linear regime, is transmitted by a normal 
fluid velocity which increases with the logarithm of the 
distance from the vortex core, but this must be cut off at 
some distance R c at which nonlinear inertial effects be- 
come more important than the viscous force. The cut-off 
may alternatively be set by the spacing between vortices, 
by the finite size of the system, or by finite viscous dif- 
fusion length at nonzero frequency of vortex motion. It 
was initially surprising to us that for these conventional 
boundary conditions there should be a transverse force, 
but it exists because the coupling between normal fluid 
flow and superfluid flow at small distances from the in- 
ner boundary produces it. We find that the magnitude 
of this transverse force is proportional to 



Pi 



m [ln(i? c /r )] 2 



(1.1) 



where u n is the normal fluid velocity at distances of the 
order R Cl tq » L is the distance from the axis at which 
inner boundary conditions are applied, and L is the mean 
free path for excitations. This transverse force due to the 
normal fluid is small unless the radius of the boundary is 
comparable with the mean free path. A significant result 
that we get from our generalization of Lamb's solution 
jl5| of the Oseen equation is that at large distances from 
the vortex there is a normal fluid circulation, contrary to 
what was assumed by TAN, but this circulation is pro- 
portional to 1/ In 2 U 7 where U is the asymptotic normal 
fluid velocity, so it goes slowly to zero as the normal fluid 
velocity goes to zero. Our results do agree with TAN's 
relation between the transverse force and the circulation 
of momentum at large distances. Some details of the 
calculation are shown in Appendix |A|. 

In sec. ^ we give a brief discussion of how the inner 
boundary conditions can be modified to allow for the 
transverse force generated by flow of normal fluid past 
a bare vortex. Our study of the asymptotic flow shows 
that any force dependent on the normal fluid velocity is 
modified by a logarithmic dependence on the cut-off dis- 
tance R c . In the low velocity limit any force due to the 
normal fluid motion, such as the Iordanskii pd[ | force, 
is renormalized away by logarithmic denominators; the 
force on a bare vortex has the form given by eq. (1.1) 
in the case ro ~ L. In practice these terms should not 
be ignored, and may give a significant normal fluid con- 
tribution to the transverse force on a vortex, in addition 
to the superfluid Magnus force. The relation between 
vortex velocity and the transverse force on it does have a 
topological structure similar to the relation between volt- 
age and current in a quantum Hall device [Em, but, unlike 



the quantum Hall conductance, this relation is subject to 
large corrections due to the finite size of the system, or 
the nonzero value of the velocity, decreasing only as the 
reciprocal of the square of the logarithm of the relevant 
length parameter. 

There is a discussion of the relation of this work to the 
original work of Hall and Vinen in Appendix [b]. 

In a recent paper by Sonin ||l7|| there is a similar dis- 
cussion of the relation of the TAN results to two-fluid 
theory. Sonin argues that there must be an additional 
contribution to the Berry phase from the flow of phonons 
or rotons past a vortex, and that this force is transmitted 
to large distances by the usual hydrodynamic processes. 



II. EQUILIBRIUM SOLUTIONS 

The equations of two-fluid hydrodynamics are given 
in the book by Khalatnikov @, in eqs. 9-13, 9-14, 9- 
15, 9-16, 8-14, 8-19' and 9-11. We want these equations 
in steady state conditions, when all time derivatives are 
zero. The equation of conservation of matter is then 



V • (p s v + p n u) = , 



(2.1) 



where we have written v for the superfluid velocity, and 
u for the normal fluid velocity. The momentum balance 
equation (Navier-Stokes equation) is 

VkipsViVk + PnUiUk) + ViP = V fc [77(VfeUi + VjUfc)] 



Ctp- C2)v • u + Ciu • Vp 



(2.2) 



where r\ is the first viscosity, and the £s are coefficients 
of second viscosity in the two-flui d m odel. We have set 
V • j equal to zero here, from eq. (2.1). The equation for 
superfluid dynamics integrates to give 

P-P0 + \v 2 = (C 4 - pCs)V • u - Csu • Vp . (2.3) 



The equation for thermal balance is 

V.(Su-«^)=f , (2.4) 

where S is the entropy density, n is the thermal conduc- 
tivity, and R is the rate of heat generation by dissipative 
processes, a term which is quadratic in the normal fluid 
velocity and temperature gradient, and will be set equal 
to zero in the linear approximation. Finally, there is a 
thermodynamic relation 



dp = SdT + pdp + ^/9 ti d(v - u) 2 



(2.5) 



The equilibrium vortex, if the pinning center has cylin- 
drical symmetry, is described very simply by the solution 
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u = , v 



dp 



V , P = P0~ ~ o~0 ! 



T = Tn 



This is c omp atible with the local thermodynamic rela- 
tion, eq. (2.5). 



III. LINEARIZED STEADY STATE EQUATIONS 



Now we want to linearize the equations fl2.1|)-(pr5|) 
about the vortex solution (|]^). We will take full ac- 
count of first order terms in Sv, u, Sp, ST, Sp, but will 
neglect the consequent variations of p s ,p n , S, and of the 
transport coefficients. This is an additional simplifying 
assumption, and does not affect the asymptotic form of 
the solutions, since such variations all come into the equa- 
tions accompanied by an additional factor of r~ 2 . We 
consider only motion in the xy plane, so all equations 
will be two-dimensional. The equation of conservation of 
matter becomes 

p„V • u = -p s V 2 w , (3.1) 

where the superfluid velocity has been written as 

v = V(T^ + w), where W = h9/m . (3.2) 

It is convenient now to write the normal fluid velocity in 
the form 



U{ ^ij^ 7X 



Pn 



(3.3) 



where is the alternating tensor in two dimensions. The 
linearized Navier-Stokes equation can then be written 

vJpsiVjW^jW + Sp) + p s (V l W)^J 2 w 



= 7 7ey -V J V 2 x-^(^-CiP + C 2 )V l (V 2 u;) . (3.4) 

Pn <J 

This can conveniently be split into two parts by taking 
the curl and divergence of the equation. The curl gives 



(3.5) 



which is an equation relating the viscous force on the 
normal fluid to the momentum flow in the superfluid. 
The divergence gives 

V 2 Sp + 2 Ps V l V J [{V t W)V ] w} 



-(^-Cip + C 2 )(v 4 w ) = o 

Pn « 



(3.6) 



The superfluid equation gives 



Ps 



6(1 + (ViW)ViW = -— (C4 - p( 3 )V 2 w 

Pn 



(2.6) The heat flow equation gives 



^l W 2 w+^W 2 ST = 0. 

Pn T 

The linearized thermodynamic relation gives 
Sp = SST + pSp + p(ViW)ViW 

- Pntijty iW)V jX ■ 



(3.7) 



(3.8) 



(3.9) 



Substitution of eqs. (3J)), (3J3), (3/7), in eq. ( p.q ) gives 
2p a V i V j [{V i W)V j w] - p, i V 2 [e lJ (V l VK)V 3 x] 



+— hv- pCi 

Pn 3 



S 2 T Pi 



pC4 + C 2 +p 2 C 3 )(V 4 u;) 



KPn 



V 2 w = 



(3.10) 



The last two terms on the left side of this equation de- 
scribe counterflow between normal and superfluid com- 
ponents, while the first two terms give coupling of this 
counterflow to the vortex motion. From now on we will 
ignore the contributions of second viscosity, which greatly 
simplifies the analysis. Only for £2 do we have experimen- 
tal values (from the observed attenuation of first sound) 
f§, although we know from Khalatnikov's @ eq. (9.12) 
and the inequality (9.18) that their contribution to the 
coefficient of V 4 w is pos itive. With this simplification 
we can rewrite eqs. ( p.5[ ) and ( 3.10 ) in a dimensionless 
form. These coupled linear fourth order equations have 
eight independent solutions for each angular symmetry. 

Since temperature gradient and heat flow may be im- 
portant both in the interpretation of these equations and 
i n th e boundary condition s, w e tak e th e gradient of eq. 
~9h, and substitute eqs. (3.4) and fl3.7j) to get 



SVtT = - Ps (V l W r )V 2 w - PsViUVjWWjw] 



— i/V l V 2 w + 77e 4J V J V 2 x 

Pn & 



+ PnV i [e ife (V j W)Vjbx] . 
We introduce the length scale L defined by 



L 



3S 2 T 



(3.11) 



(3.12) 



From kinetic theory L should be of the order of the mean 
free path for the excitations. We rewrite the equations 
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in terms of the dimensionless complex variables z = [x + 
iy)/L and its conjugate complex z, using the symbols d 
and 8 to denote the partial derivatives with respect to 
these two variables. In these variables we have 



W 



h z 

- — m — 

2mi z 



(3.13) 



for the vortex potential. In these variables there is one di- 
mensionless parameter in the equations of motion, which 



3frp n 
4mrj 



(3.14) 



this is of order unity at high temperatures, but falls off at 
lower temperatures in 4 Hc in a similar manner to p n /p, 
since Tip/mrj remains of order unity. We look for solutions 
of the form 



p s w = $tzf(zz) , p n x = ^szg(zz) 



(3.15) 



since these have the symmetry which corresponds to uni- 



form flow. We use s to denote the variable zz — ? 
and then the components of velocity are given by 

L Ps (v x - iv y ) =f + sf + sf'e- 2ie , 

Lp n {u x - iuy) = g- f + s(g' - f) 

-s(g> + r)e-™. 
In these new variables eq. (|3.5|) gives 



(3.16) 



H s ^2 + 2^-fg = -za(2s 2 /"' + 7sf" + 2f) , (3.17) 
and eq. ( 3.10D gives 



d 2 „ d l w d 2 „ d v 



ds 



ds 2 



ds' 



-2ia{sf" + /') + ia(2s 2 g"' + 5sg") 



(3.18) 



It is help ful to consi der t he case a — 0, where the two 
equations (3.17) and ( 3 . 18) ) decouple. In this case eq. 
( |3.17 ) has four solutions, namely 



goes, goclns, g constant , g oc 1/s 



(3.19) 



The first of these gives a solution of the Navier-Stokes 
equation with flow which is quadratic at large distances, 
the second is the well-known solution that is logarithmic 
at large distances and gives rise to a constant viscous 
force on any surface surrounding the axis, the third gives 
uniform flow, and the fourth dipolar flow falling off like 
1/r 2 at large distances. Equation ( 3.18| ), in the limit 
a — > 0, has solutions 



/ constant , / oc 1/s 



/exs- 1 / 2 /^) , f<xs- x l 2 K x {Js), 



(3.20) 



where I±, K\ are modified Bessel functions. The first two 
of these, combined with the corresponding solutions of 



eq. ( 3.19 ), give uniform and dipolar superfluid flow, while 
the third and fourth give counterflow of the normal and 
superfluid components, growing and decreasing exponen- 
tially with distance from the axis. 

Solutions of the actual equations ( 3.17 ) and ( ft.lS ) have 
the same character at large distances, but some of them 
are modified at small distances from the axis. The four 
solutions that are unchanged are /' = 0, g' = 0, g' con- 
stant, and the solution / = l/s,g= — 1/s. Combinations 
of the first two give uniform superfluid or normal fluid 
flow, the third gives quadratically growing normal fluid 
velocity, and is not of interest to us, while the fourth gives 
dipolar superfluid flow. The real part of the solution 



/ = PsVq , g = p s v , 



(3.21) 



gives uniform superfluid flow, no normal fluid flow. The 
temperature is a constant T everywhere, and the other 
thermodynamic variables are given by 

Sp = p s vo sin 6* , Sp — v Q sin (9 . (3.22) 

mr mr 

This is an exact solution of the equations as we have 
approximated them, but would be perturbed by terms 
from the inhomogeneity in the equilibrium value of p s . 
These perturbations lead to corrections of w which are 
of order 1/r, and to corrections of the pressure which are 
of order 1/r 3 . The net force per unit length on a large 
cylindrical volume of fluid of radius R due to pressure is 



f 2 \ Aa nTl - 

/ dprdo = p s z x Vo , 

Jo rn 



(3.23) 



while the momentum flowing into the cylinder per unit 
length is 

—p s 6v -fd6 = -—p s zxv . (3.24) 
m to 

This gives the usual expression for the Magnus force on 
the vortex pinning center as the sum of these two contri- 
butions, equal to 



F M = —p s z x v 

777 



(3.25) 



This derivation is exactly the same as the derivation for 
the Magnus force in a classical compressible fluid. 
Uniform normal fluid flow is given by 



/ = , g = p n u 



(3.26) 



This gives uniform pressure, no momentum flow, and no 
viscous forces. The dipolar superfluid flow is obtained 
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by taking the derivative of the vortex solution with re- 
spect to the position of the vortex core, which is, in our 
notation, 



h{s)=iai K f (s,t)F (t)dt 



f(s) oc i , g(s) = -f(s) 



(3.27) 



Various methods are available to find the remaining 
four solutions of the equation, or, rather, the three solu- 
tions that do not grow exponentially at large distances. 
We have solved the equations numerically, we have stud- 
ied the power series solutions in positive or negative pow- 
ers of s, and we have looked at the systematic expansion 
of the solutions in increasing powers of a using a stan- 
dard Green function technique for solving the inhomoge- 
neous ordinary differential equations that arise. Since a 
is indeed a small parameter for much of the temperature 
range of superfluid 4 He, and the expansion in powers of a 
is more transparent than the other methods, we will out- 
line this method, and mention some of the results from 
the other methods. 

We start with a zero order approximation that is a 
linear superposition of the solutions in the limit a = 0, 
excluding those that increase faster than logarithmically 
at large distances, and excluding the uniform superfluid 
flow, which we know about already, so that we have 



9o(s) = ci +c 2 hc 3 ln(s/a ) , 

s 



fo(s)=d 2 — + d 3 ^=K 1 (^) . 
s V s 



(3.28) 



Subst itution of this in the right hand sides of eqs. (3.17), 
( |3.18 ) gives two functions iaGo,iaFo which have the 
form 

Go = -ds^(VS) , 



F = -2(ca + d 2 )^ - c 3 - + ds-^K^s) , (3.29) 

S S l\j s 

where aL = tq is the distance at which short distance 
boundary conditions will be imposed. This leads to the 
inhomogeneous equations 

/ d 2 d ,2 



Particular integrals of these inhomogeneous equations 
can be obtained by a standard Green function technique, 
in the form 



gi(s) = ia K g (s,t)G (t)dt 



(3.31) 



and solutions d x +c' 2 a 2 /s, d 2 a 2 / s + d' z Ki{yfs)/ t/s, should 
be added to the particular integral so that the boundary 
conditions at s = a 2 are satisfied. We keep the coeffi- 
cient of ln(s/a 2 ), the dominant term at large distances, 
constant during the iterative procedure. The Green func- 
tions Kf,K g can be written as 



Kf(s,t) 



= / 1/s - (2/VS*)iTi (V*) if t < s, 

1/t - (2/v^/l (yffiK-L (V~t) if t > s, 



K g (s,t) = 



_ / (21ni- 



(21ns 



t/s)/6 
- s/t)/6 



if t < s, 
if t > s, 



(3.32) 



but it is actually convenient sometimes to use other forms 
which differ from these by a solution of the homogeneous 
equation. 

The process can be repeated by sub sti tutin g /i , <7i , 
etc., in the right hand sides of eqs. ( 3.17 ), ( |3. 18 ) to gen- 
erate Fi,Gi, and so to find the second order corrections 
in f,9- 

All the solutions generated in this way have series so- 
lutions in ascending powers of s, with the leading power 
no given by the solutions of the complex indicial equation 

127^(7^ - 1) + Qianl - a 2 {Anl - 1) = . (3.33) 

The solution generated by a combination of c 2 , d 2 has 
a power series solution in descending powers of s, with 
1/s as the leading term. The solution generated by C3 is 
Ins times a descending power series, and is the solution 
which leads to a net viscous force at large distances. The 
solution generated by d 3 has a singular point at infinity. 



IV. BOUNDARY CONDITIONS AND 
SOLUTIONS 

In the previous section we have shown how linearized 
solutions of the two-fluid equations in the presence of a 
straight superfluid vortex can be found. We must com- 
bine these solutions in a way that satisfies appropriate 
boundary conditions. We already know that superfluid 
flow gives rise to a Magnus force whose magnitude de- 
pends on the superfluid velocity at large distances. It 
remains to find the corresponding solution for the nor- 
mal fluid component. As with a classical fluid, we cannot 
satisfy boundary conditions simply with a flow which is 
asymptotically uniform, but we need the logarithmically 
growing term. Four relations between the other five co- 
efficients must be determined by boundary conditions on 
some inner boundary. 

If the fluid is rotating around a wire of circular cross- 
section, radius r — aL, as in the Vinen experiment [p2| , 
it is natural to take both components of normal fluid 
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velocity, and the radial component of super fluid velocity, 
to be zero at ro, or, in the notation of eq. ( [3.15 ), 



2a 2 /V) + /(a 2 ) = , g(a 2 ) = 



2a 2 g'(a 2 ) - f(a 2 ) = 



(4.1) 



An additional condition must be imposed, such as that 
the temperature is uniform on the boundary (alterna- 
tively, it might have been more realistic to tak e the heat 
flux to be zero at the boundary). From eq. (3.11) the 
condition for constant temperature at r — aL is 



Psh , 2 1 
( V w + 



maL 



3p n aL 



+ (??v 2 



-^jj^d e )d r x = , 
na z L z 



(4.2) 



or, in terms of / and g, and making use of eq. ( |4.1| ), 

3(2aV" + 7a 2 g" + 2g') = 4(1 - ia)(a 2 f" + 2/') . (4.3) 

Only the terms for which V 2 x and V 2 w are nonzero con- 
tribute to this. 

Our main concern in this work is to what extent the 
normal fluid flow past a vortex can generate a transverse 
force. The direction of the viscous force is determined 
by the phase of the coefficient of the term in g which is 
logarithmic in s for large s, while the direction of flow is 
determined also by the constant term. The existence of 
a transverse force depends on the complexity of the solu- 
tions of eqs. ( [3.17 ) and (3.18), imposed either by the form 
of the equations or by the boundary conditions. Numeri- 
cal solution shows that in this case (no tangential normal 
fluid flow or temperature gradient at the boundary) there 
is a transverse force due to normal fluid flow, and that 
this transverse force is proportional to p n and more or less 
independent of 77. This transverse force is small when the 
radius of the inner boundary is significantly greater than 
L, but rises as a approaches unity. 

In Appendix |A] some of the details of the solution of the 
equations with these boundary conditions to first order 
in a and lowest nonvanishing order valla are shown. To 
zero order in a and for a » 1, eq. (A3) gives 



ci 



[-! + £]*. 



(4.4) 



which shows that the solution is similar to the solution of 
the classical problem, and the coupling between normal 
and superfluid components simply reduces the effective 
radius of the cylinder by an amount close to 3L/4. 

There are imaginary corrections to this solution to first 
order in a, coming both from the combined effect of the 
Green functions in eq. (3.31), and from the boundary 
conditions. The detailed work in Appendix ^ shows that 
the most important terms for a » 1 are of order ia/a 2 . 



When the coefficient C3 of the logarithmic term in u is 
kept equal to a constant C3 in the iterative procedure, 
the coefficient of the constant term is, from eq. (A10), 



[-l + |§^]Cs-^C 3 . (4.5) 
1 2aAi(a) J a z 



In general the normal fluid velocity at large distances 
should have the form 



Lp n (u x 



.,) = C 3 [ln(s/a 2 ) + 1] - C 3 e~ 2ie + C x 



(4.6) 



Because of the logarithmic term, this does not tend to 
uniform flow at large distances, which is the basis of 
the well-known Stokes paradox for slow viscous flow past 
a cylinder. The resolution of the Stokes paradox was 
provided by Oseen in 1910 (see refs. JT3] and who 
pointed out that at sufficiently large distances the non- 
linear inertial term, p n u ■ Vu, in the equation of motion, 
overwhelms the viscous term, ?7V 2 u, however small the 
fluid velocity u is. This was discussed in the context of 
the drag on superfluid vortices by Hall and Vinen ]l2] | and 
by Vinen 13 1, and Appendix |B] shows how our discussion 
relates to this early work. In the region where these non- 
linear terms are important there is negligible coupling 
between the superfluid and normal components, so the 
classical theory can be taken over with little modifica- 
tion. 

Oseen proposed replacing the nonlinear term pu- Vu in 
the Navier-Stokes equation by pU- Vu, where U is a con- 
stant vector equal to the asymptotic velocity of the fluid. 
This gives the leading term in a small Reynolds number 
asymptotic expansion, uniform in r, and is a good ap- 
proximation in the regions where the nonlinear term is 
more important than the viscous term. The vorticity for 
the normal fluid then satisfies the equation 



p n U ■ curl u = 77V curl u 



(4.7) 



which follows from the curl of eq. (2.2) by dropping all 
terms involving the divergence of u. The general solution 
of this equation with the right asymptotic behavior, with 
U in the x direction, is 



curlu = e kx J2PnK n (kr)e m0 , 



(4.8) 



where the p n are arbitrary coefficents, K n is a modified 
Bessel function, and 



PnU 

2 V 



(4.9) 



Lamb |lq] gave an explicit solution for u, and a simple 
generalization of this solution to allow for a transverse 
component of the force is 

u x -iu y = U- Ae kx K {kr) - A e - l9 [e kx K^kr) - 1/kr] 



G 



U - Ae kL ( z+ *V 2 K (kLVzi) 



AJ-e kL ^/ 2 K 1 (kLV^) + 



A 

kLz 



(4.10) 



It can be verified that this indeed gives incompressible 
flow whose curl satisfies the Oseen equation, since we 
have 

— -^-[u x — iu v ) = curl u + i div u 

L dz y v> 



2ke kx [%AK - ^Ae 1 ")^} 



(4.11) 



The vorticity given by this expression vanishes for kr > > 
1, except in a parabolic wake bounded by 9 2 w 2/fcr; 
outside the wake the flow is a combination of dipolar and 
vortex motion. A general discussion of this method has 
been given by Proudman and Pearson |2J 



The leading term in the expansion of eq. (4.10) in pow- 
ers of k is 



iuytvU + A]n(kL\z\/2) + Aj - 



A Az 
~2 ~ ~2z~ 



(4.12) 



where 7 is Euler's constant. In the range L << r « 1/k 
this should match the form given in eq. (4.6), and this 
requires 



A = 2C 3 /Lp r , 



(4.13) 



and 



Lp n U = C 3 [2 ln(2/fcLo) - 2 7 ] + C 3 + C 3 + d . (4.14) 

The complex ratio C\/C 3 is dete rmine d by the solution 
of the differential equations ( 3.17 ) and ( |3.18 ), and by the 
boundary conditions, so the phase s of C \ , C 3 are deter- 
mined by the imaginary part of eq. (4.14), and the ampli- 
tudes are related to the asymptotic normal fluid velocity 
U by the real part. In the case that we have studied in 



detail this is given in eq. (4.5) 



Substitution of eq. (1.6) in eq. ( |2.2j ) gives the change 
of pressure and the viscous stress tensor in the region 
L «r « 1/k as 



Sp - 



4t) 
rLp n 



4:7] 

rLp r 



(9?C 3 cos (9-3C 3 sin 6) 



{rf-99)(?R.C 3 cos9-%C 3 sm9). (4.15) 



Pressure and viscous stress each contribute half the net 
force, so that the total force per unit length due to normal 
fluid flow is 



F„ = M(HC 3 £-9fC 3 0) , 

Lp n 



(4.16) 



in agreement with the results for an ordinary fl uid qu oted 
i n hy drodynamics text books Jl4|]. Equations ( 4.10 ) and 
(4.13) show that the normal fluid has circulation 

4tt3C 3 



kLp r , 



(4.17) 



in the region r >> 1/k. Comparison of this with eqs. 



(4.16) and (4.9) shows that the transverse component of 
the force is given by 



F1 



PnK n U . 



(4.18) 



This is in full agreement with the equations derived in 
TAN [Q, but not with the discussion, since it was as- 
sumed that the normal fluid circulation is zero. In fact, 
as we shall show, K n goes to zero as U goes to zero, but 
only as fast as 1/| lnC/| 2 



The real parts of eqs. ( 4.14 ) and ( 4.16 ) combined with 
the approximation ( |4.4| ) gives 

F„ • x 87777 



U 



(4.19) 



2 \n(Ar]/p n Ua) ~ 27 + 1 + 3/2a 

which is the standard result from classical hydrodynam- 
ics, apart from a small correction of order 1/ a in the de- 
nominator. Combination of the imaginary parts of eqs. 



(4.5) and (4.14) gives 



3C, 



4a /a 2 



-3?C 3 



(4.20) 



21n(477/p n [/a) - 1 + 3/2a 

so that, after substitution of the value of a from eq. 
( |3. 14 ) , the transverse force is given by 

3p n (h/m) 



(4.21) 



o?[\n{4T 1 /p n Ua)- 1 /2f 

which agrees with the result written down in cq. fll.l| ). 
As the cut-off R c = 4r]/p n Ua goes to infinity both these 
components go to zero, the transverse force faster than 
the longitudinal force. For any finite value of the cut-off 
these corrections of order l/lni? c and l/(lni? c ) 2 will be 
needed^ Our nu meric al results agree with these expres- 
sions fl4.19D and ( |4.2l| ). 

This calculation shows that even a solid core can give 
rise to a transverse force which is proportional to the nor- 
mal fluid density, and its ratio to the normal fluid velocity 
close to the core can be quite large if the core radius is 
comparable with L (a of order unity), however, its ratio 
to the asymptotic value of the normal fluid velocity is 
reduced by a factor proportional to l/(lni? c ) 2 . 

As was observed by Hall and Vinen, the actual cut- 
off may be reduced to lower values, if a nonzero angular 
velocity produces a finite spacing of vortices, producing 
a situation analogous to the array of cylindrical obstacles 
considered for a viscous fluid by Sewell , or if there is 
a nonzero frequency v of oscillation of the normal fluid 
relative to the vortex, in which case a cutoff of magnitude 



JL 

PnV 



(4.22) 



is appropriate. 
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V. FREE VORTICES AND THE IORDANSKII 
FORCE 

For a free vortex core, or one bound by a core whose 
size is much less than the mean free path of excitations, 
then we should take account of the effect of the super- 
fluid flow on noninteracting excitations, in accordance 
with the discussions of Lifshitz and Pitaevskii pof , Ior- 
danskii |ll| , and others . According to these works 

the flow of phonons or rotons past a stationary vortex 
produces a transverse force equal to 



F t = z x u c , 

m 



(5.1) 



where u c is the normal fluid velocity at a distance of the 
order of a mean free path from the vortex core. In the 
phonon dominated regime there is a longitudinal force 



_ Pn hk B T 15ttC(5) 
lp m mc 2 4£(4) c 



(5.2) 



where c is sound velocity and £ denotes the Riemann £- 
function. In the roton dominated region the longitudinal 
force is 



m rnksT 



2L 

]n(y/2tik B T—)u c , (5.3) 



where fi is roton mass and po is the magnitude of roton 
momentum. For helium in the phonon-dominated region, 
the longitudinal force due to normal fluid flow is about 
40% of the transverse force, while in the roton dominated 
region the longitudinal force is considerably larger than 
the transverse force. 

A detailed calculation of what happens in this case 
is not easy to undertake, since we have no satisfactory 
theory of the transition region between the collisionless 
region of the order of a mean free path in radius and the 
hydrodynamic region beyond a few mean free paths from 
the vortex core. The best we have been able to do is 
to join one region abruptly to the other. Qualitatively 
the same thing will happen that happens for a solid core, 
which is that the growth of the logarithmic term in the 
normal fluid velocity will lead to an increasing alignment 
of the direction of normal fluid flow with the direction 
of the force at larger distances from the core, until the 
Oseen or other cut-off radius R c is reached. Indeed, with 
a solid core whose radius r is comparable with the mean 
free path L, two-fluid hydrodynamics, although it is not 
strictly applicable, predicts a transverse force comparable 
with the Iordanskii force. The phonon-dominated regime 
is a little different, because the transverse force is initially 
somewhat larger than the longitudinal force, so, even af- 
ter the logarithmic growth of the normal fluid comes into 
effect, there is still a significant angle between the direc- 
tion of the force and the direction of flow. 



VI. DISCUSSION 

In this paper we have investigated how to combine the 
TAN expression relating the transverse force on a vor- 
tex to the circulation of momentum round the vortex to 
the usual results based on a combination of two-fluid hy- 
drodynamics with the Lifshitz-Pitaevskii-Iordanskii ex- 
pressions for the force on a free vortex in the presence 
of superfluid and normal flow. We have found that for 
a single stationary vortex in an infinite superfluid the 
TAN relation between transverse force and circulation of 
momentum is correct — this is not surprising, since the 
relation depends on little more than conservation of mo- 
mentum and the independence of normal and superfluid 
flows at large distances from the vortex. The assumption 
that is wrong in rcf. [I] is that there is no normal fluid 
circulation, which is only literally true when the normal 
fluid is at rest. Application of the methods of Oseen and 
Lamb to the two-fluid situation shows that normal fluid 
flow with speed U round a vortex induces a circulation 
proportional to p n hU / ' m[hi(4:irr)U / p n a)] 2 . 

For an isolated vortex this factor proportional to the 
inverse square of a logarithm will be rather small, and 
so the transverse force due to normal fluid flow will be 
small compared with the transverse force due to super- 
fluid flow. There is a possible modification of this con- 
clusion when the normal fluid component is primarily 
phonons, but then the normal fluid density is very small 
relative to the superfluid density, so the transverse force 
due to the normal fluid flow is still going to be small. 

Similar conclusions have been reached in the recent 
paper by Sonin p7||. 
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APPENDIX A: FIRST ORDER CORRECTIONS 

We work out the first order corrections to the flow, tak- 
ing the specific case of the boundary conditions in which 
there is no normal fluid flow and constant te mper ature 
at r — aL, where a is large, as given in eqs. (Rjj) and 



flyD 



To zero order the functions g and / are given by 



8 



of s of 1 
g = ci+c 2 hc 3 ln— , / = d 2 \- d 3 —=Ki(*/s) 



(Al) 



The four boundary conditions give 



1 



d 2 = d 3 K[(a) , c 2 + ci = , — 2c 2 + 2c 3 = d 2 + d 3 -Ki(a) , 3c 3 = -(1 - ia)d 3 aKi(a) 

a 

These equations have solutions 

3K (a) 



ci = -1- 



2aK 1 (a)(l - ia) 



C3 



-c 2 , d 2 



3K[(a) 



aKi{a)(l — ia) 



^■c 3 , d 3 



aKi(a){l — ia) 



c 3 



Substitution of eq. (3.29) in eqs. (3.31) and (3.32) gives 

+2 „2 



9i{s) = d 3 -j [2\n--—^]K[(Vt)dt = d 3 - 



t + s 

£572 



(A2) 



(A3) 



(A4) 



where three integrations by parts have been performed, and 



/i(s) = ia(c 2 + d 2 ) 



+-nr^ K ^)h{a) 



OO ey 2 



2a 2 (t - s) 
si 5 



dt 



Aa 2 



iac 3 



— dt 

st 



s l/2£5/2 

s 2 



2 S 



L/2^3/ 



K 1 (^)h{Vt)dt + 
Kx{V~s)h^/i) dt 



ia 2 



s l/2 t 5/2 

OO o 



,1/2^/2 



hi^K^V^dt 



h{y/a)Ki{y/i) dt 



+ ^72^2 K i(^) 7 i( a )J +*'«*{- _/ i^ K 'i(St)dt + J ■^[K 1 (^)I 1 (Vt)-I 1 (V^)K 1 (Vi)]K! 2 {Vi)dt} 



To leading orders in a, s, these give 

4 25 



ffi(s) « d 3 iaJ^( 



,3s 5 / 4 6s 7 / 4 
Now we substitute the expressions 



/i(s) w ia (c 2 + d 2 W - c 3 - ln(s/a 2 ) - d 3 



7T 1 

2T7 4 ' 



5 i (s) + ci + 4^ , A( S ) + 4y + 4^^i(V5) 



(A5) 



(A6) 



(A7) 



in the boundary conditions ( }4.1| ) and (4.3) to get 

c'x = [-§(2aV + 7aV + 2a 2 g[) + 2{a\f'{ + 2a 2 /()] + a 2 (f[ g[) + h 9i 



4-gi, d , 3 K (a) = 2c' 2 + 2a 2 (fl-g , 1 ) + 2f 1 , d' 2 = d' 3 K[(a) + 2a 2 f[ + f- 



(A8) 



The coefficient ci is the most important one, since, with ci,c 3 it dominates the behavior for large s. In the ex- 
pression for c' x it can be seen that the only terms giving contributions of order ia/a 2 are the terms proportional to 
a 3 g'{ , a 2 g[, a 3 /{' , a 2 f[, /i, while the term proportional to a 5 g'i gives a contribution of order ia/a which cancels with 
a similar term in c±. There is no contribution to this order from the terms in fi proportional to d 2 or d 3 . With the 
substitution of eqs. (A3) and (A6) in this we get 



ci , K (a) . . 1 , HF 2 _ a r 3K a (a) 4 , 



Combined with eq. (A3) this gives 



ci + cj _ r 3K (a) , _ ia 

c 3 ^ L 2aK 1 {ay a 2 



(A9) 



(A10) 



The last term in eq. (A10) gives the imaginary part of the coefficient, which leads to a velocity at large distances 
which is not in the same direction as the force. 







APPENDIX B: RELATION TO THE 
HALL— VINEN THEORY 

As is clear from earlier parts of this paper, some of the 
problems that we have addressed were foreseen by Hall 
and Vinen Q (HV) and by Vinen |13|, and it may help 
the reader if we examine the extent to which the results 
derived in this earlier work accord with those derived 
here. The earlier work referred to a vortex with a core of 
microscopic size (much smaller than the excitation mean 
free path); furthermore, as we have noted, it did not 
include any transverse force due to flow of the normal 
fluid, such a force having been ruled out, incorrectly, by 
a supposed symmetry of the scattering of excitations by 
the vortex. However, it is easy to modify the earlier work 
by the addition of a transverse force, as we now show. 

In the presence of a transverse component of the force, 
equation (2) of HV must be generalized to 



/,, = Dvm - D'v R± , f ± = Dv R± + D'v m 



(Bl) 



The symbols || and _L refer to directions parallel and per- 
pendicular to the asymptotic normal fluid velocity, U (i.e. 
/|| = F„ • x and f± = F„ • y in the notation of section 4 
), and Vn is the velocity of the normal fluid close to the 
(stationary) vortex. Equation (4) of HV still holds and 
gives the relationship between U and v# 



h 

E 



U - v m = — , -v R ± = — 



E 



(B2) 



where E is given by equation (5) of HV for the case of 
finite frequencies. The parameter L in E is roughly the 
mean free path of the excitations constituting the normal 
fluid, as in earlier parts of this paper. In the case of 
zero frequency, with which this paper is concerned, the 
par ameter A in equation (5) of HV must be replaced by 
jyl p n U/2r], so that the logarithmic term in equation (5) 
of HV becomes lni? e , where R e is the Reynolds number 
Lp n U/4ri. When U is sufficiently small the magnitude 
of this (negative) logarithmic term is much larger than 
unity, and then E can be written in the approximate form 



E 



In R P 



From these equations we can derive the relations be- 
tween the components o f the force on th e vortex and U , 
analogous to equations ( 4.19 ) and ( 4.21 ) in section 4. 



h_ E 



(D + E)E 2 



E) 2 + D' 2 



u 



E 2 D' 



(D + E) 2 + D' 2 



(B4) 



For sufficiently small values of U we must have D ^> E. 
It follows that the leading term in /y /U is proportional to 



[In R e ] 1 , while the lead ing te rm in f±/U is proportional 
to [lni? e ]- 2 . Equations ( ftljj| ) and < \%.2l\ ) in section 4 ex- 
hibit the same feature. According to HV the feature has 
its origin in the form of equations (B2), which give the re- 
sponse of the normal fluid to the line force f . The cut-off 
in the solution of the linearized Navier-Stokes equation 
that is required at large distances (section 1) forces this 
response to decrease logarithmically as U decreases and 
the corresponding cut-off distance increases. 



[2] 
[3] 

[1 

[5 
[6 

[7 

[S 
[9 

[10 
[11 

[12 

[13 
[14 

[is; 

[16 



(B3) [17 



[18 

[19 

[20 

[21 

[22 
[23 

[24] 



D. J. Thouless, P. Ao and Q. Niu, Phys. Rev. Lett. 76, 
3758 (1996). 

C. Wexler, Phys. Rev. Lett. 79, 1321 (1997). 

M. R. Geller, C. Wexler and D. J. Thouless, Phys. Rev. 

B 57, R8119 (1998). 

J. M. Tang and D. J. Thouless, Phys. Rev. B 58, 14179 
(1998). 

G. E. Volovik, Phys. Rev. Lett. 



77. 



4687 (1996). 
Rev. Lett. 80, 



4356 



H. E. Hall and J. R. Hook, Phys. 
(1998). 

E. B. Sonin, Phys. Rev. Lett. 81, 4276 (1998). 
E. B. Sonin, Phys. Rev. B 55, 485 (1997). 
R. J. Donnelly, Quantized vortices in helium II (Cam- 
bridge University Press, 1991). 

E. M. Lifshitz and L. P. Pitaevkii, Zh. Eksp. Teor. Fiz. 
33, 535 (1957) [Soviet Phys. JETP 6, 418 (1958)]. 
S. V. Iordanskii, Ann. Phys. (NY) 29, 335 (1964); Zh. 
Eksp. Teor. Fiz. 49, 225-36 (1965) [Soviet Phys. JETP 
22, 160 (1965)]. 

H. E. Hall and W. F. Vinen, Proc. Roy. Soc. A238, 215 

(1956) . 

W. F. Vinen, Proc. Roy. Soc. 242, 493 (1957). 

G. K. Batchelor, Introduction to Fluid Dynamics (Cam- 
bridge University Press, 1967), pp. 240-6. 

H. Lamb, Hydrodynamics, 6th Edition (reprinted by 
Dover, New York, 1945), pp. 615-7. 

D. J. Thouless, in Topological Aspects of Low Dimen- 
sional Systems, ed. A. Comtet, T. Jolicoeur, S. Ouvry 
and F. David (EDP Sciences and Springer), p. 790. 

E. B. Sonin, to be published in Proceedings of the Isaac 
Newton Institute Workshop on Quantized Vortex Dynam- 
ics and Superfluid Turbulence, August 2000. 

I. M. Khalatnikov, An Introduction to the Theory of Su- 
perfluidity, translated by P. C. Hohenberg (W. A. Ben- 
jamin Inc., New York, 1965). 

J. Wilks, The Properties of Liquid and Solid Helium 
(Clarendon Press, Oxford, 1967), Chapter 8. 
E. B. Sonin, Zh. Eksp. Teor. Fiz. 69, 921 (1975) [Sov. 
Phys. JETP 42, 469 (1976)]. 

M. Nielsen and P. Hedegard, Phys. Rev. B 51, 7679-99 
(1995). 

W. F. Vinen, Proc. Roy. Soc. 260, 218 (1961). 

I. Proudman and J. R. A. Pearson, J. Fluid Mech. 2, 237 

(1957) . 

C. J. T. Sewell, Phil. Trans. Roy. Soc. A210, 239 (1910). 



10 



